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, Abstract. We report on the method developed by Zibetti, Chariot & Rix (2009) to construct 
\^ resolved stellar mass maps of galaxies from optical and NIR imaging. Accurate pixel-by-pixel 
colour information (specifically g — i and i — H) is converted into stellar mass-to-light ratios 
with typical accuracy of 30%, based on median likelihoods derived from a Monte Carlo library 
of 50,000 stellar population synthesis models that include dust and updated TP-AGB phase 
I I , prescriptions. Hence, surface mass densities are computed. In a pilot study, we analyze 9 galaxies 
- spanning a broad range of morphologies. Among the main results, we find that: i) galaxies appear 
' much smoother in stellar mass maps than at any optical or NIR wavelength; ii) total stellar mass 
• , estimates based on unresolved photometry are biased low with respect to the integral of resolved 
' stellar mass maps, by up to 40%, due to dust obscured regions being under-represented in global 
colours; iii) within a galaxy, on local scales colours correlate with surface stellar mass density; 
^ iv) the slope and tightness of this correlation reflect/depend on the morphology of the galaxy. 

Keywords, galaxies: general - galaxies: stellar content - galaxies: structure - galaxies: funda- 
^ ' mental parameters - galaxies: photometry - techniques: photometric 
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The key role of stellar mass to determine (or predict) the physical properties of present 
day galaxies has been established by a number of works since the last decade (e.g. 
Gavazzi et al. 1996t Scodeggio et al. 20021 : iKaufFmann et al.l 120031 ). Indi cation has been 



provided that stel lar mass density might be even more fundamental (jBell fc de Jong 



200nt IKaufFmann et al...2003) . All these works, however, deal with global estimates of the 



^\ I stellar mass-to-light ratio, which is assumed to be uniform throughout a galaxy, contrary 
r to what we should expect based on well known spatial variations of stellar population and 
^ ' dust properties. Hence, resolving the distribution of stellar mass is crucial to properly 
measure total stellar mass and mass density and to start investigating the structure and 
dynamics of galaxies in an unbiased way. Having access to the mass distribution also al- 
^ . lows to address questions like: Is there any relation between local stellar mass density and 
" " the local SED and physical properties within a galaxy, similarly to the relations observed 
globally? If so, what does it tell us about the internal mechanisms of galaxy evolution? 

With these goals and quest ions in mind, we have devel oped a new method to build 
stellar mass maps of galaxies ( Zibetti. Chariot Sz Rixl 2009L ZCR09 hereafter), which we 



review in this contribution (Sec. 2). We also present preliminary results on the colour- 
stellar mass density relation within galaxies for a small sample of different morphological 
types (Sec. 3). 
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2. From multi-band imaging to stellar mass 

As described in ZCR09, our goal is to provide a computationally fast and observa- 
tionally cheap method to obtain stellar mass maps. We express the surface stellar mass 
density at any given position (a, S) of a galaxy as: 



(2.1) 



where Ea(q!, 5) is the surface brightness in the chosen 'luminance' band of effective wave- 
length A and T\(a,6) is the corresponding effective M/L ratio. In turn, we want to 
express as a function of colour indexes. Despite the diversity of stellar population 
parameters and dust properties, if the luminance band and the colour indexes are ap- 
propr iately chosen, the sca tter of as a function of colour(s) can be reasonably small 
(e.g. iBell fc de Jong 120011 ). We choose the H band (largely equivalent to J or K) as 
luminance band and we use g — i and i — H colour indexes to predict M/L^. This choice 
is supported by the fact that M/L variations in NIR are minimal with respect to shorter 
wavelengths, while these two colours are the best combinat ion in terms of sen s itivity 
and wavelength leverage. Based upon the 2007 version of the lBruzual fc Chariot! (|2003f l 
code, which i mplements updated prescriptions for the TP-AGB stellar evolutionary phase 
according to lMarigo et al.l (j2008r ). we build a Monte Carlo library of 50,000 stellar popu- 
lation synthesis (SPS) models with a variety of star f ormation histor i es, bo th continuous 
and bursty, metallicities and dust attenuations a la ICharlot fc Fal] (|200df). In this way 
we mean to cover as uniformly as possible the parameter space occupied by individual 
regions in galaxies on scales of ~ 100 pc, from old, metal-rich bulges/spheroids, to young 
spiral arms with different degrees of dust absorption. Then we bin the models in cells of 
0.05 X 0.05 mag^ in g ~ i and i — H and for each cell we consider the median M/L//. 
This look-up table is then used to assign the median-likelihood M/L to each pixel in a 
galaxy, depending on its colours. Our approach is therefore bayesian rather than frequen- 
tist and, as such, depends to some extent on the chosen prior distribution of models: a 
different choice of SFHs in particular can affect the M/L of the bluest stella r popula- 
tions (see ZCR09 for a detailed discussion). We confirm (e.g. iMara ston 2005 ) that the 
treatment of the TP-AGB phase in the SPS models has a great systematic impact on 
the resulting colours and M/L in NIR, but this is mainly limited to young/blue stellar 
populations. Apart from systematic effects, we find that the confidence half-range for the 
M/L is « 30% over most of the colour space, with the notable exception of the stellar 
populations with ages < 1 Gyr, where the stellar evolution is very rapid and results in 
larger ranges, up to a factor 2. The use of two colours is crucial to obtain such accurate 
estimates: by neglecting the dependence on i — H one would easily over/under-estimate 
the M/L by a factor 2 to 30 

To test our stellar mass mapping meth od we select a sample of 9 nearby galaxies, 
which are also part of the SINGS survey (jKennicutt et al. I l2003h and for which optical 



SPSS imaging and medium-deep NIR im aging (either from GOLD Mine. iGavazzi et al 



2003L or UKIDSS. [Lawrence et al.ll2007l) are available. The sample spans all range of 
morphologies, from ellipticals to late spirals. As fluxes and colours in individual pixels 
have to be accurate at few percent level in order for our method to work, it is necessary 
to pre-process the images to ensure that sufficiently high S/N can be reached also in 
the regions with surface brightness as low as 24-25 magg arcsec"^. This is done using 



f In ZCR09 we show that using i-band and the g — i colour as M/L predictor gives results 
that are in very good agreement with those based on g — i — H in most cases, except in stellar 
populations dominated by a very young burst or with very large extinction. 
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a newly developed adaptive smoothing code ADAPTSMOOTH (Zibetti in prep., ZCR09), 
which preserves the maximum spatial resolution compatible with the requested S/N. 

Figured] shows the maps of effective M^/Lj^jf/ for the 9 galaxies of the sample. Earlier 
type galaxies have more uniform and higher M/L, on average. The young stellar popula- 
tions that characterize spiral arms result in lower M/L ratios, which in turn decrease the 
spiral arm contrast. The resulting stellar mass maps reveal that the structure of spiral 
galaxies is significantly smoother than it appears at any optical or NIR wavelength. The 
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Figure 1. The maps of effective stellar mass to NIR light ratio of the 9 galaxies, sorted by 
morphological type, from Elliptical to Sc. Darker grey represents higher M/L, as given by the 
scale aside. 

colour gradients, especially visible in spiral galaxies, result in M/L gradients, which in 
turn affect the structural parameters and make them differ from the ones derived from 
light distribution. In particular, the effective radius of spirals is smaller when the stellar 
mass distribution is used. 

In Fig. [T] one can see lanes of much higher M/L (especially NGC4569 and NGC4536), 
which correspond to dust lanes. These regions can hide a significant amount of stellar 
mass, although the emerging light hardly affects the global luminosity and colours of the 
galaxy. As we show in ZCR09, this causes stellar masses obtained from global photometry 
to be biased low with respect to what one gets from integrating resolved stellar mass 
maps: the latter may miss up to 40% of the total stellar mass of a galaxy if dust obscured 
regions are very extended (as in NGC4569 and NGC4536). 



3. Trends of colours with surface stellar mass density 

Based on the analysis of the previous section, we study the correlation between colours 
(namely the optical g — i) and the surface stellar mass density within each galaxy. The 
distribution of pixels as a function of these two quantities are shown in Fig. [21 where the 
grey scale denotes the density of pixels. As a general result, we see that the colour pos- 
itively correlates with surface stellar mass density: galaxies are redder in higher density 
regions. Although this holds for all galaxies, the slope and the dispersion of the correlation 
varies a lot along the morphological sequence. For a typical Elliptical galaxy the relation 
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Figure 2. The distribution of pixels in the g — i vs surface stellar mass density plane for the 9 
galaxies: in grey scale the log of normalized number of pixels per cell. The horizontal line marks 
the global colour of the galaxy. 

is very flat and very tight (consistent with a scatter due to photometric errors only). 
As we move to later types the colour-mass relation steepens, with low-density regions 
becoming increasingly bluer, while the highest density regions have roughly constant red 
colours. We can interpret this as a sign that star formation prefers low density and that 
the colour-morphology relation is mainly set by the relative weight of the younger stellar 
populations in the lower density regions with respect to the red high-density 'cores'. In 
addition, we see that the scatter around the mean relation increases going to later types, 
which have a more disomogeneous distribution of physical properties at given surface 
mass density. The presence of dust also increases the scatter, especially at high surface 
mass density. These preliminary results based on a small sample of only 9 galaxies will be 
put on a much more solid ground, both from the statistical and physical point of view, 
by the forthcoming analysis of a larger sample of galaxies for which multiwavelength 
observations, from UV to radio, are available. 
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